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Abstract 

Static color charges in Yang-Mills theory are considered in the Schrodinger pic¬ 
ture. Stationary states containing color sources, interquark potential and confinement 
criterion are discussed within the framework based on integration over gauge transfor¬ 
mations which projects the vacuum wave functional on the space of physical, gauge- 
invariant states. 


1 Introduction 

Usually the potential of interaction between static quark and antiquark is extracted from 
Wilson loop expectation value. A Wilson loop is essentially non-local operator, since it 
describes the process developing in time. On the other hand, static charges are point-like 
and it might seem that in the Schrodinger picture, when stationary states are considered, 
they can be described by some local operators. Of course, such operators can not be local 
in terms of the original gluon fields, but it is reasonable to look for auxiliary variables, in 
which color charge creation operators can be local. The best candidates for such variables 
are pure-gauge components of the gluon fields. 

* Based on talks given at 10th International Seminar on High Energy Physics ”Quarks’98”, Suzdal, Russia, 
May 18-24, 1998 and at Workshop on ’’Instantons and Monopoles in the QCD Vacuum”, Copenhagen, 
Denmark, June 22-27, 1998. 
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Longitudinal gluons are unphysical in Yang-Mills theory, but, once color charges are con¬ 
sidered, longitudinal degrees of freedom do not decouple completely, since the electric held of 
a static charge is longitudinal and longitudinal gluons are responsible for color electric forces. 
In this respect, it is useful to introduce gauge degrees of freedom as independent variables 
from the outset. Such kind of variables naturally appears if the gauge invariance of physical 
states is imposed by averaging over gauge transformations. An approach to the ground state 
in Yang-Mills theory based on this procedure was advocated in Refs. i a. where it was 
proposed to find an approximate vacuum wave functional from variational principle applied 
to the simplest, Guassian variational ansatz. Averaging over gauge transformations allows 
to make trial wave functional exactly gauge invariant. It also allows to construct easily 
the states describing static color sources and to formulate a relatively simple criterion of 
confinement in the Hamiltonian framework [3j. £1]. 


2 Vacuum sector 


The ground state in gluodynamics satisfies Schrodinger equation 

H ym ^ = E v ac T, (2.1) 

where -Hym is Yang-Mills Hamiltonian (in A 0 = 0 gauge): 



When the wave function is a functional of the gauge potentials, the electric fields act as 
variational derivatives: 



(2.3) 


We consider SU(N ) gauge group with generators T A obeying conventional normalization 
conditions tiT A T B = *-5 AB /2 and sometimes use matrix notations for gauge potentials: 
At = A a T a . 

The Hamiltonian (|2.2|) commutes with operators DiE A (x), and, apart from the Schrb- 
dinger equation, the physical states are subject to the Gauss’ law constraint: 


D,E a T = 0. 


(2.4) 


The Gauss’ law represents in the infinitesimal form the gauge invariance of the physical 
states: 

T[M n ] = V[A], A? = (yA t + X - a (2.5) 

To be more presize, the infinitesimal form of the invariance condition allows for an additional 
factor in eq. ( [2.5|) related to topologically non-trivial gauge transformations 0 . Absence of 
this factor corresponds to the case of zero theta-parameter, which is implied hereby. A 
non-zero theta is discussed in Ref. [lj. 


2 



There are many ways to solve the Gauss’ law constraint - to use gauge-invariant variables 
[ 0 , to fix the residual gauge freedom [[7|, or to enforce gauge invariance projecting the wave 
function on the physical subspace § 0- 0, 0- 1- The projection operator has an elegant 
path-integral representation as an averaging over all gauge transformations: 


'[/[.4] = I [DU] e -'1 J P 


( 2 . 6 ) 


This representation is well known in Yang-Mills theory at finite temperature [Kj and was 


used to find an approximate ground state of the Yang-Mills theory by variational method 
[0, |2], In a sense, eq. (|2.6| ) is a general solution of the Gauss’ law, since any gauge-invariant 
functional can be represented in this form by an appropriate choice of the ’’bare” state 
exp(— S[A]). 


3 Static charges 

3.1 Coupling to heavy fermions 

To introduce static charges, we couple Yang-Mills fields to infinitely heavy fermions. By 
infinitely heavy we infer fermions without kinetic term in the Hamiltonian: 


H = fd 3 x( i E*E* + i F$F* + M#) ■ (3-1) 

The fermion operators obey anticommutation relations 

{CW, $}fb{yj) = S a pS a b 8(x -y), (3.2) 

where a,b and a,f3 are color and spinor indices, respectively. The fermion fields transform 
in the fundamental representation of SU(N ): 

4> u = (3.3) 

The ground state of the Hamiltonian ( |3.1| ) can be easily constructed, because H obeys 
simple commutation relations with fermion operators: 

[HMx)\ = ±Mi> a {x ), |+’ ° a Zl A 2 (3-4) 

= ±Mipl(x), |^’ • (3.5) 


These relations are written in the basis in which 7 0 is diagonal. It is clear that ^ 1 , 2 (^), ^3 ±(x) 
are annihilation operators and the fermion vacuum is defined by equations 

Vh,2(z)|0) = 0, 

4A x )\°) = °- ( 3 - 6 ) 
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The fermion vacuum is gauge invariant by itself, and the ground state of the Hamiltonian 
(|3.1|) is described by the wave function 


T vac = *[A]\0) = I [DU] e- 5 ^] |0), (3.7) 

where T[A] is the vacuum wave functional in gluodynamics - the lowest-energy solution of 
the Schrodinger equation ( |2.1| ). 

Operators ^^(x), ^l 2 ( x ) create excited fermion states which contain static color charges. 
These states are not gauge invariant and, thus, it is necessary to project them on the phys¬ 
ical subspace. Again, projection can be realized as averaging over gauge transformations. 
Therefore, physical states are obtained by acting of gauge transformed operators on the 
fermion vacuum with subsequent integration over the gauge group. This procedure can be 
considered as a kind of dressing of the bare fermions by their static color-electric fields. The 
dressed operators are gauge-invariant - the gauge indices are replaced by global color ones 
[[|, which we mark by prime: 

= y'kWCW. 

( 3 ' 8 ) 

where a = 3,4 and (3 = 1,2 (below spinor indices are omitted for brevity). The simplest 
example of a color-singlet state containing static charges is the one with meson quantum 
numbers: 

= I [DU] e- s[AU] if>lF{x)i/} Ua '(y)\0) = ^ a b [A;x,y] ipi(x)i/j b (y)\0), (3.9) 

where the gluonic part of the wave function, 

np ;*,!/] = JlDU]e- s ^U\,(x)U^' b (y) (3.10) 

transforms in the representation N ® N of the gauge group in order to compensate gauge 
transformations of the fermions. 


3.2 Abelian theory 

The above construction, in fact, generalizes the dressing of electron operators proposed by 


Dirac 11 


( :\y ) = e~ ieV ^My), 

U 


(* )f ^ - - ieV{x) jpUx), 


, where 


= e 

1 


V(x) = Jtl [ y —yjr. - I;,). 


(3.11) 


(3.12) 


Dressed operators (|3.11|) possess a number of important properties. They are gauge-invariant 
and create eigenstates of the free electro-magnetic Hamiltonian. They also play an important 
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role in QED eliminating IR divergencies related to emission of soft photons |[T2| . |~Hf|| . Different 
non-Abelian generalizations of the operators (|3.11|) were proposed |12], [14 . 


It is not difficult to reproduce the dressing (|3.11 ) from the representation based on 
averaging over gauge transformations. In Abelian theory U = e ie G Af = A* + and 
the action S'[A] is quadratic; for simplicity we choose the diagonal quadratic form. Then 


the integration over gauge transformations in (|3.10|) is Gaussian and, for example, the state 
describing two charges of opposite sign has the form: 


V[A;x,y] = l[dtp\exp 

= C e ieV ^ e ~ ieV( <A exp 


— 2 K S- diip) 


e ieip(x) e —ie<p(y) 


^KAf 


(3.13) 


where C is a field-independent constant and Af denotes transversal part of the gauge po¬ 
tentials: 


d 2 


(3.14) 


As follows from conformal symmetry, the coefficient function K is equal to \p\ in the mo¬ 


mentum space. Then the last factor in (|3.13|) is nothing but the vacuum wave functional 
for free electro-magnetic field. The first two factors reproduce Dirac dressing operators. For 
a generic charged state the operator ?/> u will be replaced by after averaging over gauge 
transformations due to the Gaussian nature of the integration over the Abelian gauge group. 


4 Interquark potential 

The energy of the state ( |3.9|) after obvious subtractions determines qq interaction potential: 


(T A /|Tm) 


— 2M + E vac -\rV(x — y). 


(4.1) 


Matrix elements of gauge-invariant states, like those entering eq. (|4 . Ip , are proportional to 
the volume of the gauge group. The representation (|2.6| ) allows to get rid of this infinite 
factor easily |l| . For example, the norm of the vacuum state is 

(T|T) = J[DU][DU'][dA] e -S[A^]-S[^'] = const j[DU][dA] e -^ e ]-s[^]. (4.2) 

This expression can be viewed as a partition function of a three-dimensional statistical 
system. Matrix elements of the Hamiltonian are determined by a linear response of this 
system on perturbation by the operator 


R[A\ = J d 3 x 


-1 y^4]_ _ 1 3SJ41 3S(A| 1 ' 

2 SAfSAf 2 SAA SAf 4 « «, 


(4.3) 


which, essentially, is the Hamiltonian. So, defining the partition function 


Z = 


, R[A u ]+± R[A] 


(4.4) 
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we express, for example, the vacuum energy in gluodynamics as a derivative of the free 
energy of the statistical system ([4.4Q with respect to A: 


WHymW .8 
v “ 3A 


(4.5) 


A=0 


This quantity is proportional to the volume and contains UV divergent contribution from 


zero-point oscillations. Regularized energy density can be related to gluon condensate [|15 


E = 


/o F A F A ^ V o') Vol + UV divergent terms, 
16a, ' ' ' 


(4.6) 


where /3(a s ) is the /3-function. 

The quark-antiquark potential can be represented in the form 


d 

V{x-y)= — In 


tr U(x)U\y) 


(4.7) 


A=0 


where averaging is dehned by eq. (14). The confining behavior of the qq potential is rather 
natural from the point of view of this representation. Really, the correlation function in 
eq. ([4.71) must decrease at infinity. So, its logarithm, the derivative of which determines 
the potential, always increases. This does not mean that the potential is always rising at 
infinity. If the correlator of gauge rotations increases as a power of distance with fixed 
exponent, the potential decreases. However, if a mass gap is generated and the correlator 
falls exponentially, the potential grows linearly: 


tr U{x)U\y)) oc 


r = \x — y\ 


oo, 


V(r) = or + ..., 

with the string tension determined by the derivative of the mass with respect to A: 

dm 

dX A=0 


a = — 


(4.8) 

(4.9) 


(4.10) 


Hence, the confinement arises due to generation of a mass gap in the averaging over gauge 
transformations. More presizely, the confining potential is generated if there is a nonzero 


linear response of the mass gap on the perturbation by the Hamiltonian in eq. fl4.4|) . 


4.1 Short-distance behavior 

At short distances the potential can be calculated perturbatively. The free-gluon (zeroth 
order perturbative) wave functional is Gaussian - the action S is quadratic: 

S[A}= l -A?KA*+.... (4.11) 

We do not specify the form of the kernel K, since it drops from the final answer. The 
matrices of gauge transformations are expanded as U = e 9$ = 1 + g& A T A + ..., and 
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AV = Ai + 8$ + .... The gauge potentials can be integrated out in this approximation by 
the change of variables: Ai = Ai — di&/ 2, and the partition function (|4.4|) takes the form 


Z = const 


exp 


-\ d$ A (a + 1a 2 ) 


1 


(4.12) 


For the correlator which defines the interquark potential we get: 

^ tr U(x)U^(y)^ = N + —/$ A (x)$ A (y) — - & A (x)$ A (x) 

2 \ 2 

= N+ S AAZZL(D(x-y)-D(0))+O(g l ), 

where D is the propagator of the held <3>: D ~ 1 = —d 2 (K + \K 2 /2)/2. The equation ( f4.7|) 
now yields the Coulomb potential with correct coefficient: 


(4.13) 


V(x-y) = 


g 2 (N 2 - 1) / 1 


2N \-d 2 
g 2 (N 2 - 1) 


[x 


y) 


8N7TT 


+ self-energy. 


4.2 Large-distance behavior 


The large-distance asymptotics of the potential is determined by infrared structure of the 
vacuum wave functional, which, generally speaking, is unknown. For illustrative purposes 
we present here a calculation of the string tension under particular assumptions about the 
dependence of the wave functional on long-wavelength fields ||. 

The starting point is derivative expansion, which probably is justified at large distances. 
The main assumption is that the derivative expansion of the kernel K in eq. (|4.11|) starts 
from the constant term (in momentum representation): K(p ) = M + 0(p 2 ). This assumption 
is compatible with variational estimates of Ref. |l| ; however, arguments in favor of M = 0 
also exist [|TS|. If we accept that M ^ 0, cubic and higher terms in S [4], as well as the 
magnetic part of the Hamiltonian, are irrelevant being of higher order in derivatives. Thus, 


the action in (|4.4|) takes the form of a local expression quadratic in gauge fields: 


Z = 


J [dA} [DU] exp 


-- (m + -m 2 
2 l 2 


d 3 x (AfA 


% + A 


U A_/\U A 


(4.14) 


The integration over gauge potentials can be eliminated by the change of variables Ai = 
Ai — diUU^/2g. The remaining path integral over gauge transformations is the partition 
function of three-dimensional principal chiral held with local action: 

Z = f [dA] [DU] exp 


■ 2 _ (m + - M 2 
2 g 2 \ 2 


d 3 x trdjWdjU 


(4.15) 
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Principal chiral field is expected to produce a mass gap [1IJ • The simplest way to see it is 
to consider U and U' as independent complex matrices and to impose the unitarity condition 
introducing a Lagrange multiplyer. If the Lagrange multiplyer acquires non-zero vacuum 
expectation value, the field U has massive propagator. Following Ref. |l]], we use mean field 
approximation to determine the mass gap. The gap equation follows from unitarity condition 

UW) = 1: 


d 3 p 1 
(27r) 3 p 2 + m 2 


2 g 2 N 


M + ^M 2 


(4.16) 


Differentiating this equation in A we get, according to ( |4. 10|) , 


dm 7i M 2 M 2 

~d\ = g 2 N = Aa s N' 


(4.17) 


This relation is obtained within the effective low-energy theory. Hence, the coupling here 
is not the running one; it is fixed on the typical energy scale: ot s = a s (M). It is worth 
mentioning that the corrections to the zeroth order of derivative expansion and to the mean 
field approximation are parametrically not suppressed and can only be numerically small [|lj . 


5 Discussion 


The averaging over gauge transformations appears to be very convenient for consideration of 
static color charges in the Schrodinger picture. It allows to formulate rather simple criterion 
of confinement. But originally integration over the gauge group was proposed in the context 
of the variational approach |l| . In principle, variational calculations can be done analytically 
for the simplest Gaussian ansatz projected on the space of gauge-invariant states. The main 
problem in this approach is correct UV renormalization. If UV divergencies are not removed 
by standard counterterms, a variational approximation looses sense. The UV structure of 
the vacuum wave functional was discussed from different point of views in Refs. |T7 . 


Another interesting question concerns topologically non-trivial gauge transformations. 
Kinematical consequences of large gauge transformations are well known || . The dynamics 
of topologically non-trivial configurations of the field U was considered recently |U| for the 
Gaussian variational ansatz averaged over the gauge group [|18|] . The skyrmion in this model 
was interpreted by the authors of Ref. []I8| as a Hamiltonian counterpart of the instanton. 
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